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A brief overview is presented of a new Caltech/Cornell research program tliat is exploring the 
nonlinear dynamics of curved spacetime in binary black hole collisions and mergers, and of an initial 
project in this program aimed at elucidating the flow of linear momentum in black-hole binaries 
(BBHs). The "gauge-dependence" (arbitrariness) in the localization of linear momentum in BBHs 
is discussed, along with the hope that the qualitative behavior of linear momentum will be gauge- 
independent. Harmonic coordinates are suggested as a possibly preferred foundation for fixing 
the gauge associated with linear momentum. For a BBH or other compact binary, the Landau- 
Lifshitz formalism is used to define the momenta of the binary's individual bodies in terms of 
integrals over the bodies' surfaces or interiors, and define the momentum of the gravitational field 
(spacetime curvature) outside the bodies as a volume integral over the field's momentum density. 
These definitions will be used in subsequent papers that explore the internal nonlinear dynamics of 
BBHs via numerical relativity. This formalism is then used, in the 1.5PN approximation, to explore 
momentum flow between a binary's bodies and its gravitational field during the binary's orbital 
inspiral. Special attention is paid to momentum flow and conservation associated with synchronous 
spin-induced bobbing of the black holes, in the so-called "extreme-kick configuration" (where two 
identical black holes have their spins lying in their orbital plane and antialigned) . 



I. INTRODUCTION: MOTIVATION AND 
OVERVIEW 

A. Motivation 

Since the spectacular breakthrough by Pretorius [l| 
in spring 2005, numerical relativists have been success- 
fully simulating the inspiral, merger and ringdown of bi- 
nary black holes (BBHs). Much effort is now going into 
extracting physical and astrophysical information from 
these simulations. 

Almost all of this effort takes an "S-matrix" view- 
point: For chosen initial conditions (the two holes' initial 
masses, vectorial spins and orbital elements), what is the 
final emitted gravitational waveform and what is the final 
hole's mass, vectorial spin, and kick velocity? 

Equally interesting, it seems to us, are the things these 
simulations can teach us about the nonlinear dynamics 
of curved spacetime. This paper is the first in a new 
research program by the Caltcch/Cornell relativity and 
numerical-relativity research groups, aimed at exploring 
nonlinear spacetime dynamics in BBHs. 

B. Momentum Flow in Black-Hole Binaries 

Several sets of analytical tools already exist for explor- 
ing fully nonlinear space-time dynamics, for example dy- 
namical horizons [2] and quasi-local energy/momentum 
and angular momentum Q . One of our goals is to develop 
additional analytical and quasi-analytical tools and use 
them to extract physical insights from numerical simula- 
tions. Our initial focus in this direction is on the distribu- 
tion and flow of linear momentum in strongly nonlinearly 
curved spacetimes — with linear momentum defined via 




FIG. 1: Extreme-kick configuration for a black-hole binary: 
Identical holes, A and B with masses ?n — M/2 move in a 
circular orbit with their spin angular momenta Sa and Sb 
antialigned and lying in the orbital plane. 



psudotensors, which arise from viewing general relativity 
as a nonlinear field theory in a flat auxiliary spacetime ^ 
This paper is the first in a series that will deal with this 
subject. 

An instructive example is the extreme-kick configation 
in which two identical, spinning black holes are initially in 
a (quasi-)circular orbit, with oppositely directed spins ly- 
ing in the orbital plane (Fig- HI)- As Campanelli, Lousto, 
Zlocliower and Merritt (henceforth CLZM) discov- 

ered and Gonzalez et. al. Q helped flesh out, of all initial 
configurations, this one has the largest kick speed for the 
final black hole, ^ and it also exhibits intriguing orbital 
motions: 

During the pre-merger inspiral, as the holes circle each 
other, they bob up and down (in the z direction of Fig.[T]), 
sinusoidally and synchronously. After merger the com- 



^ Although as Chen, Nester and Tung have shown, that various for- 
mulations of pseudotensors could also be motivated from quasi- 
local points of view Q|. 

^ For binaries in non-circular orbits, larger kick velocities have 
been observed by Healy et al. Q. 
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FIG. 2: Bobbing and kick of binary black holes in the 
extreme-kick configuration of Fig. [T] as simulated by Cam- 
panelli, Lousto, Zlochower and Merritt (CLZM) Plotted 
vertically (as a function of time horizontally) is the identi- 
cal height z of the two black holes, and then transitioning 
through merger (presumably at t/M ~ 170), the height of 
the merged hole, above the initial orbital plane. This height 
versus time is shown for six different initial configurations, 
each leading to a different orbital phase at merger. In all six 
configurations, the initial holes' spins are half the maximum 
allowed, a/m = 0.5. The height z and time t are those of the 
"punctures" that represent the holes' centers in the CLZM 
computations, as defined in their computational coordinate 
system, which becomes Lorentz at large radii. These z and 
t axe measured in units of the system's total (ADM) mass 
M ~ 2m. 



bined hole gets kicked up or down with a final speed that 
depends on the orbital phase at merger (relative to the 
spin directions). This bobbing then kick, as deduced by 
CLZM from numerical simulations, is graphed quantita- 
tively in Fig. [21 

Momentum conservation dictates that, when the holes 
are moving upward together with momentum -I- , 
there must be some equal and opposite downward mo- 
mentum in their gravitational field (in the curved space- 
time surrounding them); and when the holes are mov- 
ing downward, there must be an equal and opposite 
upward field momentum. How is this field momentum 
distributed, and what are the details of the momentum 
flow between field and holes? To what extent is the fi- 
nal kick of the merged hole (e.g. in configuration SP2 of 
Fig. an incrtial continuation of the holes' immediate 
pre-merger bobbing, and correspondingly to what extent 
is the burst of downward gravitational-wave momentum 
that accompanies the kick caused by near-zone, bobbing 
field momentum continuing "inertially" on downward af- 
ter merger? And to what extent are other momentum- 
flow processes responsible for the motion shown in Fig. [2P 
These are the kinds of questions we would like to answer 
by an in-depth study of momentum flow in BBHs. 



C. Gauge-Dependence of Momentum Flow; 
Landau-Lifshitz Formalism 

The momentum distribution and fiow in a relativistic 
binary are tricky concepts, because momentum conserva- 
tion arises from, and requires, translation invariance of 
spacetime. Spacctime is translation invariant when flat, 
but not, in general, when curved. Two key exceptions 
are: (i) Spacetime is locally translation invariant in the 
vicinity of any event, and this leads to the local law of 
4- momentum conservation T°'f^.p = (where T"'^ is the 
total stress-energy tensor of all nongravitational particles 
and flelds). (ii) Around any isolated system, e.g. a BBH, 
spacetime can be idealized as asymptotically translation 
invariant, and this leads to the definition and conser- 
vation law for the system's total momentum (e.g., the 
binary's final kick momentum is equal and opposite to 
the momentum carried off by gravitational waves). How- 
ever, inside the binary the curvature of spacetime pre- 
vents one from defining a globally conserved momentum 
density and fiux in any generally covariant way. 

Nevertheless, we are quite hopeful that momentum 
flow can be developed into a powerful tool for physical 
intuition into BBHs, and into the nonlinear dynamical 
behavior of curved spacetime that is generated by col- 
lisions of spinning black holes. To do so, however, will 
require living with the fact that the momentum distribu- 
tion and momentum flux inside a binary cannot be gen- 
erally covariant, i.e. they must be, in some sense, gauge- 
dependent. 

There, in fact, is a long and successful history of physi- 
cists' building up physical intuition with the aid of gauge- 
dependent concepts; and it is that history that gives us 
hope. For example, in Maxwell's flat-spacetime electro- 
dynamics, the vector potential satisfies the wave equation 
only if one first imposes Lorenz gauge; and our physical 
intuition about electromagnetic waves relies, to a con- 
siderable extent, on Lorcnz-gauge considerations. Simi- 
larly, in developing post-Newtonian ephemeridcs for the 
solar system, celestial mechanicians have chosen a spe- 
cific gauge in which to work, and their intuition about 
relativistic effects in the solar system relics to a great 
extent on that gauge's gauge-dependent constructs. The 
choice of gauge was, to some extent, arbitrary; but once 
the choice was made, intuition could start being built. As 
a third example, in black-hole perturbation theory rela- 
tivists have built up physical intuition based on Regge- 
Wheeler gauge, and based on the Teukolsky equation, 
each of which are gauge-dependent constructs. 

The density, flux and conservation of linear momen- 
tum, in the curved spacetime of a black-hole binary, must 
rely explicitly or implicitly on a mapping of the binary's 
curved spacetime onto an auxiliary, translation-invariant 
flat spacetime. This reliance is spelled out explicitly in 
a reformulation of general relativity as a nonlinear field 
theory in flat spacetime presented in Landau and Lif- 
shitz's Classical Theory of Fields Q. (See also Chap. 20 
of MTW [l^l and a more elegant, covariant formulation of 
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the formalism developed by Babak and Grishchuk [Hj.) 
In the original Landau-Lifshitz formulation, one chooses 
any asymptotically Lorentz coordinates that one wishes, 
one maps onto an auxiliary flat spacetime by asserting 
that these chosen ("preferred") coordinates are globally 
Lorentz in the auxiliary spacetime (so in them the aux- 
iliary metric has components diag[— 1, 1, 1, 1]), and one 
then reformulates the Einstein equations as a nonlinear 
field theory in the space of that flat, auxiliary metric. 
The result is a total stress-energy tensor 

r"'' = (-,9)(r"'3 + t^{^), (1.1) 

where g is the determinant of the covariant components 
of the physical metric, T"^ is the nongravitational stress- 
energy tensor, and is the "Landau-Lifshitz pseu- 
dotensor" . By virtue of the translation invariance of the 
auxiliary spacetime, this r"'^ has vanishing coordinate di- 
vergence t""^ ,13 = in the chosen "preferred" coordinates. 
Equivalently, this t°'^ has vanishing covariant divergence 
t"^\I3 = with respect to the auxiliary flat metric. The 
components t^^ then represent the density and t^^ the 
flux of a conserved linear momentum. 

We envision each numerical relativity group choosing 
the coordinates used in its simulations to be the "pre- 
ferred" coordinates of this mapping to flat spacetime, 
resulting in each group's adopting a different "gauge" . If 
we are lucky, this will lead to momentum distributions 
and flows in different groups' simulations that are qual- 
itatively and semi-quantitatively similar. If that is not 
the case, then we advocate that the community adopt, as 
a communally-agreed-upon "preferred" coordinate sys- 
tem (and thence gauge) , Harmonic coordinates — though 
even then it might be necessary to face up to the fact 
that Harmonic coordinates are not uniquely defined un- 
til one gives appropriate initial conditions. We envision 
joint numerical and quasi-analytical explorations, over 
the coming months, that lead simultaneously to a choice 
or choices of "preferred coordinates" for the mapping to 
flat spacetime, and physical insights into the flow of mo- 
mentum in BBHs. 

This paper represents a first small step in this direc- 
tion: To ensure that we understand quite clearly what is 
going on, we shall focus in this paper on a binary's pre- 
merger bobbing, and we shall study it and its momentum 
flow using the post-Newtonian approximation to general 
relativity in Harmonic coordinates. Subsequent papers 
in this series will use the Landau-Lifshitz formalism to 
explore momentum flow in black- hole mergers. 

D. Overview of this Paper 

We begin our post-Newtonian analysis in Sec. |TT] by 
presenting our main ideas and results in the simplest in- 
teresting context: the extreme-kick configuration. 

We then, in the remainder of the paper, present a de- 
tailed post-Newtonian analysis of spin-induced momen- 
tum flow in the inspiral phase of generic compact binaries 



(BBHs, neutron-star binaries, or neutron-star / black- 
hole binaries). This detailed analysis begins in Sec. IIIII 
with a very brief summary of the Landau-Lifshitz formal- 
ism, followed in Sec. llVl bv a use of the formalism to give 
a general treatment of 4-momentum conservation for a 
fully relativistic system of compact bodies. We express 
the binary's total 4-momentum (as measured gravitation- 
ally by distant observers) as the sum of the 4-momenta 
of its two bodies (expressed as integrals over their sur- 
faces or, for stars, volume integrals over their interiors) 
and the 4-momentum of their external gravitational field 
(expressed as a volume integral over the exterior). We 
also derive expressions for the rate of change of the 4- 
momentum of each body as a surface integral of the flux 
of 4-momentum being exchanged between the body and 
the external fleld. 

In Sec.|V]we specialize to the inspiral of a generic com- 
pact binary, as analyzed in Harmonic coordinates at lead- 
ing nontrivial post-Newtonian order (1.5PN for the ef- 
fects of spin) ; and we focus on the distribution and flow 
of linear momentum (the spatial part of 4-momentum) 
induced by the bodies' spins. We begin in Sec. IV Al bv 
computing the spin-induced perturbation of the field mo- 
mentum 5t°^ in terms of the binary's masses, vectorial 
spins, and geometry; and we then integrate this density 
over the exterior of the bodies to obtain the total field- 
momentum perturbation (Jpfioid hi terms of the bodies' 
masses, spins and vectorial separation. In Sec. IV Bl we 
discuss the definition of a body's center of and 
corresponding velocity v = dxcm/dt, and we write down 
the infiuence of the bodies' spins MdSv/dt on their equa- 
tions of motion. In Sec. IV Gl and Appendix IB] we use our 
definition of center of mass to deduce an expression for 
the spin-induced perturbation of a body's momentum 5p 
in terms of its mass times velocity perturbation MSv, 
and cross products of the bodies' spins with their sepa- 
ration vector. Finally, in Sec. IVD] we verify momentum 
conservation; i.e., we verify that, as the binary evolves 
and momentum is fed back and forth between the bod- 
ies and the field, the bodies' equations of motion ensure 
that the spin-induced perturbation of the total momen- 
tum (bodies plus field) is conserved. 



II. BOBBING AND MOMENTUM FLOW IN 
THE EXTREME-KICK CONFIGURATION 

In this section we shall present an overview of our 
momentum-flow ideas and results in the context of the 
extreme-kick conflguration (Figs. [T] and [2]). 

Pretorius has offered a lovely physical explanation 
for the holes' bobbing (Fig. ^ in this conflguration: In 
Fig- El taken from his paper, we see snapshots of the holes 
at four phases in their orbital motion. In each snapshot, 
each hole's spin drags space into motion (drags inertial 
frames) in the direction depicted by gray, semi-circular 
arrows. In phase B, hole 1 drags space and thence hole 2 
into the sheet of paper (or computer screen); and hole 2 
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drags space and thence hole 1 also inward.'^ In phase D 
each hole drags the other outward. This picture agrees 
in phasing and semi-quantitatively in amplitude with the 
bobbing observed in the simulations (Fig. (2]). 

Semi-quantitatively but not quantitatively. In addi- 
tion to frame dragging, there is a second influence of 
the holes' spins on their motions, at the same 1.5 post- 
Newtonian (1.5PN) order: a force on each body due to 
the coupling of its own spin to the Riemann curvature 
tensor produced by the other body. For the extreme- 
kick configuration (Fig. [T]) , in Harmonic coordinates the 
mass times coordinate acceleration produced by frame 
dragging (expressed as a weak perturbation, "(5" , to the 
motion of a nonspinning binary) is {cPSxA/dt^)Fr> = 
{4:/r\g)SB X vab, and that produced by spin-curvature 



easily compared with Fig. [2l 



couphng is {d 6xA/dt^)sc 



r\g)SA X vab (first and 



second lines of Eq. (4.11c) of [IJ]). Here tab is the sepa- 
ration between the two holes and vab = va — vb = 2va 
is the coordinate velocity of hole A relative to hole B. 
Since Sb = —Sa, the sum of the frame-dragging accel- 
eration and spin-curvature-coupling acceleration is 



d^Sx/ 
dt^ 



-Sa X VA 



(2.1) 



spin effects 



' AB 



Wc get the acceleration of hole B by replacing all sub- 
script B's by subscript A's. The two holes' accelerations 
are identical (synchronous bobbing) because Sb = —Sa, 
Vb = -Va- 

We can easily integrate this equation in time by not- 
ing that the spin processes much more slowly than the 
orbital motion so Sa can be approximated as constant, 
and noting that va rotates with angular velocity Q = 
•y/ M/r^Q = •\/2m/r^^, where M is the total mass and 
m is the mass of each hole. The result, after one integra- 
tion, can be written as 



mSvA ~ mSvB 



-Sa X Uab ■ 



(2.2) 



' AB 



Here Sva and 6vb are the holes' velocity perturba- 
tions (bobbing) produced by their spins, and tiab — 
(xa — XB)/rAB is the unit vector that points from hole 
B toward hole A. (As with the convention for propaga- 
tors, we regard things as "moving" from right to left; the 
vector Uab points from B to A.) One might think that ex- 
pression (|2.2p represents the holes' bobbing momentum, 
but as we shall see in Sec. lHBl it does not — for subtle but 
physically understandable reasons. The bobbing momen- 
tum of each hole is actually 2/3 times expression (|2.2p : 
see Eq. (P?T5)) below. 

We can integrate Eq. (|2.2[1 once more to obtain each 
hole's bobbing displacement (change of location) in Har- 
monic coordinates. We write the result in a form that is 



5xa Sa 

= -VA X — 

m TO^ 



(2.3) 



Because Sa remains approximately constant while va ro- 
tates uniformly in time (if we ignore radiation-reaction- 
induced inspiral), and because Sa and va both lie 
in the orbital plane, Eq. (j2.3p represents an approxi- 
mately sinusoidal bobbing orthogonal to the orbital plane 
{z direction), with (peak-to-peak) amplitude 5z/M = 
vaSa/'tii^ = \vA, where we have used the spin mag- 
nitude SaI^ — a/m = 0.5 of the CLZM simulations. 

The CLZM simulations cover only the last two orbits 
before inspiral, when the Post-Newtonian approximation 
is failing badly and the inspiral is rapid. Nevertheless, we 
can hope for rough quantitative agreement. The simula- 
tion shows a maximum bobbing amplitude 5z/M ~ 0.4, 
which agrees with our 1.5PN amplitude 5z/M = ^va if 
Va is near the speed of light, as it should be just before 
merger. Half an orbit earlier the simulation's bobbing 
amplitude is smaller by about a factor 1/1.7 ~ 1/V2.5, 
which is what our 1.5PN formula predicts if the orbital 
radius is 2.5 times larger than at maximum amplitude - 
and this agrees fairly weU with Fig. 2 of CLZM ^. 



A. Field Momentum in the Extreme-Kick 
Configuration 

In Harmonic gauge at leading post-Newtonian order, 
the Landau-Lifshitz formalism gives for the density of 
field momentum 



gxg 

47r 



(2.4) 



(Eq. (4.1a) of [il])- Here, to the accuracy we need, g is 
the Newtonian gravitational acceleration field (the gravi- 
toelectric field), H is the gravitational analog of the mag- 
netic field (the gravitomagnetic field), and Bj is the j'th 
basis vector of the flat-spacetime field theory that we are 
using. (As we shall see, and as is discussed in Ref. [lH| 
and references cited therein, the analogy with electrody- 
namics can be very powerful in building up insight into 
gravitational momentum density and flux.) 

[Side Remarks: Before applying Eq. (|2.4p to the 
extreme-kick configuration, let us build up a bit of physi- 
cal insight into it: As for a particle, so for the relativistic 
gravitational field, we can regard the ratio of the momen- 
tum density to the mass-energy density, t^-' /t^" as a field 
velocity 



Afield 



(2.5) 



^ This is very similar to the way that two fluid vortices (e.g. an 
aerofoil's starting and stopping vortex pair) drive each other into 
motion; see, e.g. [Tsll . 



One can show that in Harmonic coordinates the vacuum 
field momentum density t^'^ is negative; in fact, it is 



no 
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-g g 



(2.6) 
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FIG. 3: Pretorius' physical explanation for the holes' bobbing in the extreme-kick configuration. 



at leading PN order. ^ Correspondingly, the gravitational 
field 's velocity ( as "seen" in our auxiliary flat spacetime ) 
points in the direction of +g x H and has a magnitude of 
order \H\/\g\. The direction of this field velocity is the 
same as the direction of motion of an inertial point mass 
(relative to our Harmonic coordinates) that is induced 
by a brief jomi action of g followed by H: The geodesic 
equation, in Harmonic coordinates and for low particle 
velocities v takes the Lorentz-force form 



dv 

-=g+VxH 



(2.7) 



(Eq. (2.7) of [15[) at leading order. In a very short time 
interval t, the field g acting on a particle initially at rest 
produces a velocity v ~ gt, and then H acts on this 
velocity to produce Sv = ^g x Ht^ — which points in 
the direction of the field velocity.] 

Now let us study the field momentum for the extreme- 
kick black-hole binary. For the moment we are only in- 
terested in that portion of the field momentum which is 
induced by the holes' spins, since this is the portion that 
must flow back and forth between the field and the bob- 
bing holes in order to conserve total momentum. This 
portion arises from one hole's gravitoelectric field g cou- 
pling to the spin-induced part of the the other hole's 
gravitomagnetic field 



Here 



An 



m 



rrspin 

gs X if/ 
An 



{3nA ■ SA)nA - Sa 



(2.8) 

(2.9a) 
(2.9b) 



(Eqs. (2.5) and (6.1) of [ia|), with ua the unit radial vec- 
tor pointing from the center of hole A to the field point. 



and rA the distance from the center of hole A to the 
field point. The gravitoelectric field (|2.9ap (actually the 
Newtonian gravitational acceleration) has identically the 
same form as the Coulomb electric field of a point charge, 
with the charge replaced by the hole's mass wla and the 
sign reversed. Similarly, the gravitomagnetic field (|2.9bp 
is identical to the dipolar magnetic field of a point mag- 
netic dipole, with the magnetic moment replaced by the 
opposite of twice the hole's spin, i.e., —2Sa- The fields 
for hole B arc the same as Eqs. (|2.9p . but with each sub- 
script A replaced by a i?. 

Combining Eqs. (|2.8[) and (|2.9[) we obtain for the bi- 
nary's density of field momentum (that portion which 
must flow during bobbing ^) 



[3(S^ • n^)(nA X n^) - (S^ x n^)] 
B), (2.10) 



and integrating this over the space outside the holes, we 
obtain for the total field momentum that flows during 
bobbing (the part of the field momentum that depends 
on the holes' spins)^ 



(2.11) 



This is equal and opposite to the sum of the holes' 
bobbing momenta 5pA + ^Pb, as we shall see in Sec. IH CI 
below. 

Figure 13] shows the z-componcnt (perpendicular to the 
orbital plane) of the field-momentum density (5r°^, as 
measured in the orbital plane at four different moments 
in the binary's orbital evolution. Only that part of the 
momentum which flows during bobbing [Eq. (j2.10p ] is 
pictured. Red depicts momentum density flowing out of 



This is t"'' in Harmonic coordinates (gauge) in vacuum, at lead- 
ing (Newtonian) order; see, e.g. the first term in Eq. (4.4a) of 
. This Newtonian gravitational energy density is gauge de- 
pendent; see, e.g., the discussion in Box 12.3 of [17( . 



^ There are portions of the momentum density that do not flow 
during bobbing, which will be important for our comparisons 
with numerical relativity. The full expression, therefore, is listed 
in Appendix IAI 

® Equations l|2.10| l and H2.11|l are special cases of Eqs. I|5.3ll and 
1 15.61 1 below, where the details of the integration are carried out. 
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FIG. 4; The four pictures show the ^-component of field- 
momentum density St"'^ in the orbital plane at four different 
times, a quarter orbit apart. Red represents positive momen- 
tum density (coming out of the paper), and blue, negative 
(going into the paper). Only the piece of momentum den- 
sity (5r°^ that flows during bobbing [Eq. (|2.10|l ] is depicted. 
The yellow arrows are the black holes' vectorial spins; the 
large, black arrowed circle shows the orbital path of the two 
holes. In the top-left picture, one sees the density of momen- 
tum when the black holes are at the top of their bob (max- 
imum z) and momentarily stationary [Eqs. (|2.2|l and (|2.3p ]. 
The gravitational-field momentum is zero, but the momentum 
density itself shows rich structure. A quarter orbit later, in 
the top right, the holes are moving downward (into the paper) 
at top speed. The momentum between the black holes (blue 
region) flows into the paper with them, while surrounding mo- 
mentum (red region) flows out of the paper {+z direction). 
A half orbit after the first picture, in the lower left, the holes 
are momentarily at rest at the bottom of their bob (minimum 
z), the net field momentum is zero, and the momentum dis- 
tribution is opposite that in the first picture (as one would 
expect during sinusoidal bobbing). Similarly, three quarters 
of the way through the orbit, in the lower right, the holes 
have reached their maximum upward speed, and the momen- 
tum distribution is identical to the second figure, but with 
the opposite sign. 



the paper {+z direction), and blue, into the paper. The 
yellow arrows show the holes' veetorial spins S, and the 
arrowed circle is the binary's orbital trajectory. In the 
top-left and bottom-right frames, the black holes are mo- 
mentarily stationary at the top and bottom of their bob- 
bing [cf. Eqs. (|2.2p and (|2.3p ]. Nevertheless, the momen- 
tum density has a non-trivial distribution. In the top- 
right and bottom- left frames, the black holes are mov- 
ing downward and upward, respectively, with maximum 
speed. In both cases, the field-momentum density be- 
tween the two holes flows in the same direction as the 



bobbing, whereas the momentum surrounding the binary 
is in the opposite direction and larger. This leads to net 
momentum conservation for the binary, as discussed in 
Section |VDl 

It is worth noting that the four figures, going counter- 
clockwise from the top- left, are taken a quarter period 
apart in orbital phase. The first and third differ by half 
an orbital period (as do the second and fourth); and, con- 
sequently, the momentum patterns of each pair are iden- 
tical, but signs are reversed (red exchanged with blue, 
as dictated by the symmetry of the configuration). This 
feature is responsible for the sinusoidal bobbing. 



B. The holes' momenta 

In this section we shall use a roundabout route to ex- 
plain, physically, why the momentum pA of black hole A 
is not mvA, and to derive an expression for it. 

Begin by considering, for pedagogical purposes, a 
rigidly and slowly rotating body in flat spacetime with 
rotational velocity v,:ot{x). Let p be the mass-energy 
density of the body's material in the local rest frame 
of a bit of material. Then in an inertial frame where 
the body is at rest except for its rotation (the body's 
"momentary inertial frame" ) , its mass-energy density is 
= p(l + v^ot): where ^f^ot comes from kinetic en- 
ergy and ^^rot from Lorentz contraction. We define the 
body's center-of-mass location by 

Af ... = / body's momentary rest frame , 

(2.12) 

where M = J T°^<Px is the body's mass. If the body is 
weakly gravitating, this location will be the center of the 
monopolar part of its gravitational field. 

Now let this rotating body move with a linear velocity 
V that is small compared to its rotational velocity, so 
TOO ^ + ^ „^^^)2] ^ ^ „2^^ ^ 2t;„t • v]. If we 
use this TOO compute / T^'^xd^x, we will not get the 
Xcm of Eq. (|2.12p because the term 2pt>rot ■ v will weight 
X extra heavily on the side of the body where v,-ot - v > 
and less heavily on the side where t^rot - v < 0. If we want 
to compute the correct Xcm by an integral performed in 
a frame where the body is moving, we must correct for 
this effect. The correction factor is well known ([l^ and 
Sec. ED below): 

MXcrn = J T^^xd^x -V X S , (2.13) 

where S is the body's angular momentum. Other defini- 
tions of center-of-mass are sometimes used, but they all 
differ from the locations one would identify, in the body's 
rest frame, as the mass-energy-weighted location (|2.12p 
and the center of weak monopolar gravity — i.e., they are 
less physically motivated than this one. 

Equation (j2.13p , when extended into general relativity 
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in the obvious manner, 

Ma;cm = j T°"xd^x -vx S (2.14) 

[with t"^ the total stress-energy tensor of Eq. (ll.lP l , is 
called the physical spin supplementary condition [18| : of. 
Eq. (|5.10p below, where a formal derivation is presented. 
In general relativity this condition guarantees that in 
the body's local rest frame, Xcm is at the center of the 
monopolar part of the body's (possibly strong) gravita- 
tional field, or more precisely the center of the monopolar 
part of the time-time component of the metric density 
0°*' = \/—gg'^^, which plays a major role in the Landau- 
Lifshitz field-theory-in-flat-spacetime formalism. 

The black-hole velocities va and vb used in this 
paper and in the standard Harmonic-coordinate, post- 
Newtonian equations of motion, are the (coordinate) time 
derivatives of the holes' centers of mass: 

VA^dXAcm/dt. (2.15) 

By specializing Eq. (|2.14p to body A, differentiat- 
ing with respect to time, using the conservation law 
T°*'^o + ''"'^"'0 = 0; E^J^d integrating by parts on the volume 
integral, we obtain [Eq. (|5.14p below and its derivation] 

Pa = + (^■A X SaY 

kinetic SSC 

+ / i^' - <n.A)ir"' - r'%\)dE, ,(2.16) 

JdA 

" V ' 

surface 

where 

p{= [ T^'^d^x (2.17) 

J A 

is the total 4-momentum of body A. Here is the 
acceleration of body A produced by the gravity of body B 
and m is the mass of body A. For a black hole the linear 
momentum must be defined via a surface integral rather 
than the volume integral J^r'^^d^x [Eq. (j4.1bp below], 
but Eq. (|2.16p still turns out to be true; see the paragraph 
following Eq. (|5.14p below. 

Equation (|2.16[) has a physical interpretation that is 
closely related to the one for the center-of-mass equa- 
tion (|2.14p that underlies it: Rearranged, this equation 

says mv^ ^ Pa ^ i^A x SaY - JgAi^^ - xi^A)i'^°'' - 
T°°z)^)(iEfe. The left side is the time derivative of the 
center-of-mass location, weighted by the body's mass (or 
the kinetic momentum). The first term on the right side 
is the body's total momentum, i.e. the volume integral 
of t'^^ . The second (SSC) term corrects for the fact that 
for a spinning body r^^ weights the center of mass too 
heavily on the side of the body where the rotational ve- 
locity and linear velocity are coaligned {va rot • I'A > 0) 
and too lightly on the side where they are antialigned. 



The third (surface) term corrects for a contribution to 
the momentum arising from mass flowing into and out of 
the body (mass flux t*^* — r^^t;') at different locations on 
the body's face. 

C. Momentum Conservation for Extreme-Kick 
Configuration 

We now specialize Eq. (|2.16p for the momentum of a 
body in a generic binary to the extreme-kick configura- 
tion and focus on the spin-dependent piece of the momen- 
tum that is exchanged between the bodies and the field. 
By inserting the expressions Sr^^Sj = —Qb x Ha/^t^ and 
o-A = Qsix = xa cm) hito Eq. (|2.16p . with g and H given 
by Eqs. (j2.9p . and performing the integrals, we obtain 

2 TO 

SpA = -0^2— ^A X UAB ■ (2.18) 

''"ab 

This spin-induced perturbation of the linear momentum 
of body A, when added to an equal amount for body B 
gives —^{m/r'^g)SA x uab, which is equal and opposite 
to the spin-induced perturbation of the field momentum 
(|2.1ip . Therefore, as the holes circle each other, momen- 
tum flows sinusoidally back and forth between the holes 
and the field, with no change in the total momentum; the 
total momentum is conserved during the bobbing. 

Let us examine in detail how momentum conservation 
is achieved in the presence of the bodies' bobbing. Our 
detailed analysis (above) breaks each object's momentum 
perturbation 6pa.b into three terms, the kinetic momen- 
tum mSvA.B, a term due to the SSC, and a surface inte- 
gral term (see Table |T| . The total kinetic momentum 

mSvA + m6wB = -{sla x 8,4 + as x S_b) 7^ (2.19) 

is not conserved because of the non-cancellation between 
the frame-dragging and spin-curvature coupling terms. 
The total body momentum SpA + SpB is not conserved 
cither; it sums up to 2/3 the total kinetic momentum: 

SpA + SpB = -^(a^ X + as X 85) ^ 0, (2.20) 

To achieve momentum conservation, there is a non-zero 
spin-dependent total field momentum distributed outside 
of the bodies, with 

2 

^Pfioid = g(aA X Sa + as X Ss) . (2.21) 

Note that this total external field momentum is only 
—2/3 the spin-dependent total kinetic momentum — in- 
stead of the —1 that one might have naively expected. 

It is important to notice that for each body, a canonical 
momentum can be formed by adding the SSC term to the 
kinetic momentum 

Pa canonical = mVA + ttA X Sa ■ (2.22) 
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Body 


Kinetic 


SSC 


Surface 


Total 


Frame-Dragging 


Spin-Curvature 


Pa 


-4as X Ss 


3aA X Sa 


Ha X Sa 


2 

--aA X Sa 


--aA X Sa 


Pb 


-4aA X Sa 


3as X Ss 


as X Ss 


--as X Ss 


--as X Ss 


Pficld 










2 

-(aA X Sa + as X Ss) 



TABLE I: Spin-dependent, time- varying pieces of body and field momenta at 1.5PN order, for the extreme-kick binary (circular 
orbit with spins antialigned and in the orbital plane). The body momenta are broken down into kinetic, SSC and surface terms 
and are expressed in terms of the bodies' spins Sa,b and Newtonian-order gravitational accelerations aA.s = —mnA,B/rAB- 
See Eqs. (|5.11|l and (|5.17p for a similar decomposition in a generic binary. 



The total canonical momentum pA + ps is conserved, 
because the sum of the surface terms of the bodies' mo- 
menta and the external field momentum arc equal and 
opposite. (This is so not solely for our extreme-kick 
configuration, but also for any generic binary; see Sec. 
FVl below) . This canonical momentum can be motivated 
quite simply by special relativistic kinetics, without the 
need for any knowledge of field momentum, and it is used 
in the Hamiltonian approach to post-Newtonian dynam- 
ics 

Although the introduction of canonical momentum re- 
solves the issue of momentum conservation at the level 
of two-body dynamics, it does not provide information 
about the distribution of field momentum nor the role of 
field momentum in momentum conservation. Our anal- 
ysis reveals substantial spin-dependent field-momentum 
outside of the bodies — with the same order of magni- 
tude as the total spin-dependent kinetic mometum. 

One might question the meaningfulness of a distiction 
between the bodies' (localized) momenta and the (dis- 
tributed) field momenta, because a different choice of 
gauge might move momentum between these two com- 
ponents, and possibly even move all the field momen- 
tum into the interiors of the objects. We argue that our 
choice of Harmonic gauge, the analogue of Lorcnz gauge 
in electrodynamics, is a promising tool for analyzing com- 
pact binaries, since its metric perturbations (both phys- 
ical and gauge) propagate at the speed of light, which 
will make gauge effects behave causally just as do phys- 
ical effects. For this reason, and because of the physical 
intuition that the above analysis brings, we advocate us- 
ing Harmonic gauge and its nonzero field momentum in 
analyzing compact, inspiraling binaries. 



III. THE LANDAU-LIFSHITZ FORMALISM IN 
BRIEF 

We turn, now, to a detailed analysis of momentum flow 
in generic compact binary systems. Wc begin in this 



section with a very brief review of the Landau-Lifshitz 
(LL) formulation of general relativity as a nonlinear field 
theory in flat spacctime 

This formulation starts (as discussed in Sec. HIl above) 
by introducing an (arbitrary) coordinate system in which 
the auxiliary flat metric takes the Minkowski form 77^^ = 
diag(— 1, 1, 1, 1). Gravity is described, in this formula- 
tion, by the physical metric density 

0^'' = V^g""" ■ (3.1) 

Here g is the determinant of the covariant components of 
the physical metric, and g^'^ are the contravariant com- 
ponents of the physical metric. In terms of the supcrpo- 
tential 

^M"^/3 ^ gA'^'g"/' _ gf^^'g'^P , (3.2) 

the Einstein field equations take the field-thcory-in-flat- 
spacetime form 

H^'"'"^,^f3 = IGTrr'^'' . (3.3) 

Here t''" = {-g){Tf"' + t^l) is the total effective 
stress-energy tensor introduced in Eq. indices after 

the comma denote partial derivatives (covariant deriva- 
tives with respect to the flat auxiliary metric), and the 
Landau-Lifshitz pseudotensor t'^ (actually a real tensor 
in the auxiliary flat spacetime) is given by Eq. (100.7) 
of LL or equivalently Eq. (20.22) of MTW By 
virtue of the symmetries of the supcrpotcntial (which are 
the same as those of the Ricmann tensor) , the field equa- 
tions in the form (j3.3p imply the differential conservation 
law for 4-momentum 

T^\, = , (3.4) 

which is equivalent to T^^-i, = (where the semicolon de- 
notes a covariant derivative with respect to the physical 
metric) . 

It is shown in LL and in MTW that the total 4- 
momcntum of any isolated system (as measured gravi- 
tationally in the asymptotically flat region far from the 
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system) is 

rfot = Y^£i?^"°^adS,, (3.5) 

where (Ej is the surface-area element (defined, of course, 
using the flat auxiliary metric) , and the integral is over an 
arbitrarily large closed surface S surrounding the system. 
This total 4-momentum satisfies the standard conserva- 
tion law 

= - / T«dl]j . (3.6) 

[The proof of this given in LL and MTW relies on an 
assumption that the interior of S be simply connected, 
i.e. that it not contain any black holes. However, that 
assumption is not necessary: Differentiate Eq. (|3.5p with 
respect to t, then use iJ''"°J\aO = Hf'"''\av - i?^"''^afe ■ 
The first term is — IGttt^^ by virtue of the the field equa- 
tions (j3.3p and the antisymmetry of the superpotential on 
its last two indices; and its surface integral gives the right 
side of Eq. (j3.6p . That same antisymmetry on the second 
term — i?^"'°-' ,Qfe permits us to write it as the curl of a 
3-vcctor field, whose surface integral vanishes by virtue 
of Stokes' theorem. The result is Eq. p.6p .] 

IV. 4-MOMENTUM CONSERVATION FOR A 
FULLY NONLINEAR COMPACT BINARY 




FIG. 5: The regions of space around and inside a compact 
binary system. 

is the gravitational field's 4-momentum in the surround- 
ing space. 

If cither of the bodies has a simply connected inte- 
rior (is a star rather than a black hole), then we can use 
Gauss's theorem and the Einstein field equations (|3.3p 
to convert the surface integral (j4.1bp for the body's 4- 
momentum into a volume integral over the body's inte- 
rior: 

p\= j T^^'d^X . (4.1d) 
J A 



We now apply this LL formalism to a binary system 
made of black holes and/or neutron stars; see Fig. [S] 
We denote the binary's two bodies by the letters A and 
B, and the regions of space inside them by these same 
letters, and their surfaces by dA and dB. For a black 
hole, dA could be the hole's absolute event horizon or 
its apparent horizon, whichever one wishes. For a neu- 
tron star, dA will be the star's physical surface. We 
denote by £ the region outside both bodies, but inside 
the arbitrarily large surface S where the system's total 
momentum is computed. (In later papers in this series, 
S will sometimes be the outer boundary of a numerical- 
relativity computational grid.) 

By applying Gauss's theorem to Eq. (|3.5p for the bi- 
nary's total 4-momentum and using the Einstein field 
equation (|3.3p . we obtain an expression for the binary's 
total 4-niomentum as a sum over contributions from each 
of the bodies and from the gravitational field in the re- 
gion £ outside them: 

Ptot=PA+PB+PLid- (4- la) 

Here 

p'X^^i H^"°\ad^j (4.1b) 

loi' JdA 

is the 4-momentum of body A and similarly for body B, 
and 

pLm^ j/'d?^. (4.1c) 



By an obvious extension of the argument we used to 
derive Eq. (|3.6p for the rate of change of the binary's 
total 4-momentum, we can deduce from Eq. (j4.1bp the 
corresponding equation for the rate of change of the 4- 
momentum of body A: 

^ = - / (r'^'= - r^'v\)dllu ■ (4.2) 
JdA 

Here the second term arises from the motion of the 
boundary of body A with coordinate velocity v\ = 
dx\ ^^/dt. Equation (|4.2p describes the flow of field 4- 
momcntum into and out of body A. 

We shaU use Eqs. ((iTT|) . ([^ . and speciahzcd to 

linear momentum (index /i made spatial) as foundations 
for our study of momentum flow in compact binaries. 

The actual values of the body and field 4-momenta, 
computed in the above ways, will depend on the arbitrary 
coordinate system that we chose, in which to make the 
auxiliary metric be diag[-l,l,l,l] and in which to perform 
the above computations. This is the "gauge-dependence" 
discussed above. In the remainder of this paper we shall 
choose Harmonic coordinates, so the gravitational field 
satisfies the Harmonic gauge condition 

Q'^^p = , (4.3) 

and we shall specialize the above equations to the 1.5 
post-Newtonian approximation and use them to study 
momentum flow during the inspiral phase of generic 
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compact binaries. In future papers we shall use the 
above equations, combined with numerical-relativity sim- 
ulations, to study momentum flow during the collision, 
merger, and ringdown phases of compact binaries. 



V. POST-NEWTONIAN MOMENTUM FLOW 
IN GENERIC COMPACT BINARIES 

A. Field Momentum Outside the Bodies 

As for the extreme-kick configuration, so also in gen- 
eral, the portion of the field momentum that is induced 
by the bodies' spins and that flows back and forth be- 
tween the field and the bodies, as the bodies move, is 



„ . , rrspin _ , , rrspin 



in 



47r 



(5.1) 



[Eq. (|2.8|l ]. We find it convenient to rewrite the bodies' 
gravitoelectric and gravitomagnetic fields (|2.9p as 



9'k 



rriK 



1 



2S'k ( — 



(5.2) 



where K \s A or B and where, as before, rK is the (flat- 
space) distance of the field point from the center of mass 
of body K. Inserting Eqs. (|5.2p into expression ()5.ip 
and manipulating the derivatives, we obtain the following 
expression for the field momentum density: 



(5.3) 

Here {A ^ B) means the same expression with labels A 
and B interchanged. Notice that this expression for the 
momentum density is the curl of a vector field; or, equally 
well, it can be viewed as the divergence of a tensor field. 

The total spin-induced, flowing field momentum is the 
integral of expression (|5.3p over the exterior region £ (cf. 
Fig. [5]). Using Gauss's law, that volume integral can be 
converted into the following integral over the boundary 
of £ 



Eq. (|5.4p into the form 



ZTT 



1 

dA K^A 



d'Er. 



dB \rB 



TBj .1 



[A^B). 



(5.5) 



We presume (as is required by the PN approximation) 
that the bodies' separation is large compared to their 
radii. Then on dA, we can write (1/r^) ^ = —n\/r\ 
and {\/rB),i = '^ab/'^ab' where ua is the unit vec- 
tor pointing away from the center of mass of body A, 
n^g is the unit vector pointing from the center of mass 
of body B toward the center of mass of body A, and 
Tab is the (flat-spacetime) distance between the two 
bodies' centers of mass. The first integral in Eq. (|5.5p 
then becomes uab/tb Iqa^'a/''^a^^p- ^'^^ simplicity we 
take the surface of integration to be a sphere imme- 
diately above the physical surface of body A and ig- 
nore the tiny contribution from the region between that 
sphere and the physical surface. On this sphere, we write 
dYip = r\n^dilA, where dflA is the solid angle element, 
and we then carry out the angular integral using the re- 
lation Jgj^n\n^dflA = (47r/3)(5gp. Thereby we obtain 
for the first integral in ()5.5p {4:'!T/i)SqpnAB /f's indepen- 
dent of the radius of the sphere of integration. (If 
the body is not spherical, the contribution from the tiny 
volume between our spherical integration surface and the 
physical surface will be negligible.) Evaluating the sec- 
ond integral in Eq. (|5.5p in the same way, and carrying 
out straightforward manipulations, we obtain for the ex- 
ternal field momentum 



<5pficid = 



{niBSA - mASB) X uab 



(5.6) 



For the extreme-kick configuration, which has tua = 
TUB = m and Sb = —Sa, this field momentum becomes 
expression (|2.1ip . 



Centers of Mass and Equation of Motion for the 
Binary's Compact Bodies 



27r .igs 



1 

TB 



-{A^ B) 



(5.4) 



The boundary of £ has three components: the surface S 
far from the binary on which we compute the binary's 
total momentum, and the surfaces dA and dB of bodies 
A and B. The integral over S vanishes because the in- 



tegrand is cx 1/r and the surface area is 



and S is 



arbitrarily far from the binary, r — > oo. When integrat- 
ing over the bodies' surfaces, we shall flip the direction of 
the vectorial surface element so it points out of the bodies 
(into £), thereby picking up a minus sign and bringing 



Restrict attention, temporarily, to a body that is a 
star rather than a black hole, and temporarily omit the 
subscript K that identifies which body. Then, following 
the standard procedure in special relativity (e.g. Box 5.6 
of MTW [l3|), we define the star's center-of-mass world 
line to be that set of events x^^^ satisfying the covariant 
field-theory-in-flat-spacetime relationship 

5"% = . (5.7) 

Here = / r^^d^x is the body's 4-momentum and 



— 



[(•^ ■^cm) 



r/30 



(5.8) 
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is the body's tensorial angular momentum. Here the in- 
tegrals extend over the star's interior, and because the 
star's momentum is changing, we take the time compo- 
nent of x^'jj-j to be the same as the time at which the 
integral is performed, x'^j^ = x^. (If the momentum were 
not changing, this restriction would be unnecessary; cf. 
Box 5.6 of MTW.) 

In a reference frame where the body moves with ordi- 
nary velocity = p' /p°, Eq. (|5.7p says = S'^^Vj. We 
wish to rewrite this in a more illuminating form, accu- 
rate to first order in the velocity v. At that accuracy, 
we can evaluate S'^ in the body's rest frame, obtaining 
gij ~ e^^^Sk where Sk is the body's spin angular mo- 
mentum 



Sk^ / ekiUx' -x^jT'^'^d^x . (5.9) 



Using definition ([EHD of S"° with x^jjj = x''\ our definition 
(|5.7p of the center of mass then takes the concrete form 



Thorne and Hartle] 

■niA—T— = [6nAB{SB ^ nAB ■ Vab) + 'iSs >^ VAB 

at 

-6{Sb X nAB){vAB ■ nAB)] 
-f [6nAB{SA X nAB ■ vab) + SSa x vab 

'^AB 

-iiSAX nAB){vAB ■ tiab)] ■ (5.11) 

Here 

VA = -T, , SVA = T, , VAB =Va-Vb 

dt dt 

(5.12) 

are the velocity of (the center of mass of) body A, the 
spin-induced perturbation of that velocity, and the rela- 
tive velocity of bodies A and B. The first two lines of Eq. 
(|5.1ip are due to frame dragging by the other body (body 
B)] the last two lines are a force due to the coupling of 
body A's spin to B's spacetime curvature. 

C. The Momenta of the Binary's Bodies 



majcm = / XT^^d^x -vx S . (5.10) 



Here on the left side we have replaced p'^ ~ J r'^^d'^x by 
its value in the body's rest frame, which is the the mass 
TO, since the two differ by amounts quadratic in v. 

Notice that, when computed in the body's rest frame 
so V = 0, the center of mass is mXcm = ^ XT^^d^x, hut 
when computed in any frame moving slowly with respect 
to the rest frame, this expression must be corrected by the 
term —v x S that we discussed physically in Sec MI B\ [Eg. 
\2.13\) l We asserted and used Eq. (|5.10p . the "physical 
SSC" , in our analysis of the extreme-kick configuration 
[Eq. (EH]. 

In our Harmonic coordinate system and at the 1.5PN 
order of our analysis, the dominant, time-time com- 
ponent of the Einstein field equations (|3.3p reduces to 
j^Mi'gOO _ ^gTrr^". The type of analysis carried out in 
Sec. 19.1 of MTW [l3| then reveals that in the star's rest 
frame, the monopolar part of its is centered on the lo- 
cation Xc-m') or, equivalently, when one expands the star's 
Q^^ around Xcm in its own rest frame, there is no dipolar 
1/r^ term (no mass dipole moment). This well-known 
result (e.g fl2. l2dl]) can be used as an alternative defini- 
tion of Xcm — a definition that works for black holes as 
well as for stars. 

Using this monopolar-field-centered definition of Xcm, 
Thorne and Hartle [U] have employed matched asymp- 
totic expansions (valid for black holes) to derive the equa- 
tions of motion for a system of compact bodies, e.g. a 
compact binary [their Eqs. (4.10) and (4.11)]. For a 
compact binary, the spin-induced contributions to these 
equations of motion at 1.5PN order arc [Eq. (4.11c)] of 



As in the previous subsection, we initially restrict our- 
selves to a body that is a star; then we shall generalize 
to a black hole; and we initially omit the star's label A 
or B for ease of notation. 

For a star we can derive an expression for the mo- 
mentum p' ~ J t'^^ d^x (with the integral over the star's 
interior) in terms of the star's velocity = dx^^/dt by 
differentiating the center-of-mass equation (|5.10p with re- 
spect to time. To allow for the possibility that the mass 
might change with time, we set m = J r^d'^x before 
doing the differentiation; i.e., we differentiate 

xi^ I T^'^d'^x = I x'T°°d^x -{vx Sy . (5.13) 

J A J A 

Using ft = —T°'',k and Gauss's theorem, we bring 
the left side into the form Tood^x — x{^ IdA^^^'' ~ 
T^'^v^)dYij. The last term arises from the motion of 
the surface of the star through space with velocity v. 
Manipulating the time derivative of the integral on the 
right side of Eq. (|5.13p in this same way, we bring it 
into the form /^rO^d^a; - /^^x^XrO'^ - T°°v^)dY.k = 

P'' ~ SdA^^^'^^'^ ~ T°°v'')d'Sk, where p^ is the star's mo- 
mentum. Inserting these expressions for the left side and 
the right-side integral into Eq. (|5.13p . noting that the 
star's spin angular momentum evolves (due to preces- 
sion) far more slowly than its velocity, denoting the time 
derivative of its velocity by dv/dt = a (acceleration), 
solving for p^ and restoring subscript A's, we obtain 

p>^ = mAv\+ [ (a;^ -<„^)(r"'=-r«%)dl]fc 

JdA 

+ {aA X SaY . (5.14) 

Although we have derived this equation for a star, it 
must be true also for a black hole. The reason is that 
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all the quantities that appear in it arc definable without 
any need for integrating over the body's interior, and all 
are expressible in terms of the binary's masses and spins 
and its bodies' vectorial separation, in manners that are 
insensitive to whether the bodies are stars or holes. To 
illustrate this statement, in Appendix iBl we deduce ()5.14p 
for a black hole, restricting ourselves to spin-induced por- 
tion of the momentum that is being exchanged with the 
field, Sp-'j^. 

It is this (5p^ that interests us. Because the spin has no 
influence on at the relevant order (which is (5r"*' ~ gH 
and Jr^" ^ where g and H are the gravitoclcctric and 
gravitomagnctic fields), Eq. (|5.14[) implies that 



TTlASvi 



dA 



)5r°'=dSfc + {aA X SaY 



■'cm A 



(5.15) 

The acceleration ga of body A is, at the order needed, 
just the gravitoelectric field of body B at the location 
of A, UA = — (jn 3/^^0)71 AB- Performing the surface 
integral on a sphere just above the body's physical surface 



we can write 



■'cm A 



and dTji 



rid^A- 



Inserting these into Eq. (|5.15p . we obtain 



netic terms can be read off the equation of motion (|5.1ip : 



dSvA 



dSvf 



= -{itibSa - mASs) X [vab - 3(nAs • i'ab)was] • 

By inserting ua b = (a^cmA - x cm b) I tab into the sec- 
ond term of Eq. (|5.18p and differentiating with respect to 
time, we obtain the negative of expression (|5.19p . There- 
fore, 



dSptot /dt = 



(5.20) 



i.e., as the binary's evolution drives spin-induced momen- 
tum back and forth between the bodies and the field, the 
total momentum remains conserved, as it must. 

Interestingly, during the summation of momentum 
terms, one finds that the surface terms in SpA + Sps have 
exactly cancelled the field momentum (5pfioid, leaving the 
total momentum as the sum of the bodies' kinetic term 
and their SSC term — i.e. leaving it equal to the bodies' 
total canonical momentum (see the discussion at the end 

of scc irrc]) . 



mASv\ 

kinetic 
term 



dA 



I ^11 A 



dflA 



(Sa X UabY 



' AB 



surface term 



SSC term 

(5.16) 

Here "SSC term" refers to the "spin supplementary con- 
dition" required to get the correct, physical center of 
mass; see text following Eq. (|2.14p . In the surface term, 
the field momentum density (5t°'^ is given by Eq. (|5.3p . 
The second term {A ^ B) is smaller than the first by 
M/tab and thus is negligible. Inserting the first term 
into the integral, using {\/rA),qp = {Sn'^j^n^ — Sn 



■^qp)/''~A 



and (l/rs),; 



/r^ , and J n 



bring Eq. (|5.16p into the form 
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surface term 

1 TUB 



SSC term 



Sa X Uab 



(5.17) 



VI. CONCLUSION 

In this paper we have explored the flow of momentum 
between a compact binary's bodies and their external 
gravitational field (spacctime curvature), at 1.5PN order, 
during the binary's orbital inspiral. In subsequent papers 
we shall explore momentum flow in numerical-relativity 
simulations of a binary's collision, merger and ringdown. 
We expect these studies to give useful intuitive insights 
into the internal dynamics of binary black holes and the 
nonlinear dynamics of curved spacctime. 
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D. Momentum Conservation 

The total spin-induced momentum perturbation, 
Sptot = 5pA + SpB + Spfieid [Eqs. (|5.17p and (|5^ ] is 

Sptot = mASvA+mBSvB + ^—{mBSA~mASB)^nAB ■ 

''^AB 

(5.18) 

Momentum conservation requires that the time deriva- 
tive of this Sptot vanish. The time derivative of the ki- 



APPENDIX A: THE TOTAL PN MOMENTUM 
DENSITY 

In Section HTXl Eqs. (|2?8)) and ((2?T0)) show the portion 
of the field momentum that generates bobbing. There 
are, however, additional pieces of field momentum at the 
same PN order that do not contribute to the bobbing, 
and these expressions become important for comparisons 
of the post-Newtonian analysis with numerical-relativity 
results (to be presented in future papers). There are 
three extra sources of terms. First, the gravitomagnctic 
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field has a part that depends upon the body's ve- 

locity. Using Eqs. (2.5) and (6.1) of [l^, one can see 

that 



H 



velo 



4mA (n^ X va) 



(Al) 



Second, there are terms from the coupling of the grav- 
itomagnetic field of body A with its own gravitoelectric 
field r'-^^Ej = —{qa X Ha)/{^t^), and similarly for body 
B, for both the spin and velocity pieces of H . Finally, 
there is a part due to (3f/Arg')/(47r), where J/at is the 
Newtonian potential and the dot denotes differentiation 
with respect to time (see Eq. (4.1) of [3). When one 
accounts for these additional expressions, the full field 
momentum density is written most concisely as 



spin 



clo' 



(A2a) 



where t^^-^^ and T^ly^ are the terms that depend upon the 
spins and the velocities, respectively. These terms are 
given by 



[3(Sa • nA)(nA x ns) - (S^ x hb) 



and 



j"0 



-^^(Sa xnA) + (A^B), (A2b) 



ruA f me [4(nB • va)!!^ - 4(nA • nB)vA] 



3mB(nA-VA)nB niA[{-n.A ■ ^a)^a - ^^a] 



-{A^B). 



(A2c) 



In the body of this paper we have confined attention to 
the first line of Eq. ()A2bp and its [A B), as that is the 
part of the field momentum that gets exchanged with the 
bodies during bobbing. 



APPENDIX B: MOMENTUM OF A BLACK 
HOLE 



In the text we derived expression (|5.14p for the mo- 
mentum of a body in a binary, assuming the body is a 
star so we could do volume integrals, and we then as- 
serted that this expression is also valid for black holes. 
The spin-induced portion of this expression that gets ex- 
changed with the field as the body moves is given by Eq. 
(|5.15[) , which reduces to (|5.17|) . In this appendix we shall 
sketch a derivation of Eq. (j5.17p directly from the surface- 
integral definition (j4.1bp of a black hole's momentum. 



^^'^ ^ 1^ 



6W 



aOk 



dA 



(Bl) 



To evaluate this surface integral up to desired 1.5PN- 
order accuracy turns out to require some 2.5PN fields. 



Qualitatively, this can be anticipated because the super- 
potential we use in the surface integral is sourced by the 
spin-orbit piece of field momentum, and therefore nec- 
essarily a non-leading PN term. One can see this more 
clearly by expanding SH^^'^'^^a in terms of the metric 
density and using the symmetries of the supcrpotcntial 
H (which arc the same as the Riemann tensor). In gen- 
eral, the momentum is given by 



1 

16^ 



(fl-'''0"° ■ 



ydT, 



k- 



(B2) 



dA 



In Harmonic gauge, however, q"' ,p = 0, and the spatial 
metric is flat until 2PN order while the time-space com- 
ponents are of 1.5PN order. As a result, the terms at 
lowest and next-to-lowest PN order are contained within 
two terms. 



5p\ 
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16^ 



(0' 



k- 



(B3) 



OA 



In this expression, the momentum arises from linear 
terms involving the metric density, instead of quadratic 
ones. As a result, one must keep pieces of the metric 
perturbation that are of higher PN accuracy. (Note: if 
we were to evaluate the time derivative of SpA using the 
surface integral ()4.2|) . we would not face such a delicacy; 
the integrand there is quadratic and requires only 1.5PN 
fields for its evaluation.) 

To find the momentum in terms of the standard post- 
Newtonian parameters, we resort to a standard way 
that the metric perturbations arc written in recent post- 
Newtonian literature, e.g., by Blanchet, Faye, and Ponsot 

m- 

goo = -1 + 2U-2U2 + 8A (B4a) 

5.0 = -W,-8R, (B4b) 

g^J = 6,^(1 + 2V + 2V^)+ AW,, (B4c) 

= 1 + 2V + + 2Wkk. (B4d) 

For spinning systems, we adopt the notation of Tagoshi, 
Ohashi, and Owen |22| where 0{m, n) means to order 
c™ for non-spinning terms and xc" for terms involving 
a single spin x- (Here x ~ \S\/m?' is the body's dimcn- 
sionless spin.) In this notation, terms we are interested in 
are of the order 0(3, 6), while the above post-Newtonian 
potentials have been obtained up to the following or- 
ders fH-IH: 



V 

W,k 



0(2,5), Vj 
0(4,5), Rj 



0(3,4), 
= 0(5,6). 



(B5) 







00 



In terms of these post-Newtonian potentials, U, K;, Ri, 
X and VKy, the perturbed metric density is 

1 - 4U - 2 {Wkk + 4U2) + 0(6, 7) (B6a) 

4U, -sfi?, + +0(6,7) (B6b) 

0(6,7). (B6c) 



q'' = % - 4 
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As a consequence, Eq. (|B3p is given by 
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16^ 



9A 



-Wj - 8 (^i?. 



V(,n)V.j(s) 



k 

(B7) 



where a subscript (5*) means keep only the parts of those 
potentials proportional to the spins of the bodies, and 
a subscript (m) involves pieces of the potential without 
spins (proportional to the masses of the bodies). Terms 
without a subscript have both pieces. 

Tagoshi, Ohashi and Owen express the potentials V(m) , 

Vj, Rj{s) and Wjk[s) terms of the bodies' masses, 
vectorial velocities, vectorial spins, and vectorial sepa- 
rations, and distance to the field-point location [their 
Eqs. (Ala), (Aid), (Alf), and (Alg)]. While the full 
equations are quite lengthy, the portions that generate 
momentum flow - those involving the coupling of the 
mass of one body to the spin of the other - are some- 
what simpler. For convenience, we give these portions of 
the equations below, rewritten in our notation, with the 
typos noted by G. Faye et al [2^ corrected. 
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h ejkii)A S rij^ 

'A 



(B8a) 
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AB 



3m f 



i-AB' 



^rATAB 



(B8b) 



Here, as before, m^, va, and Sa are the mass, velocity, 
and spin angular-momentum of object A; is the sep- 
aration of body A from a point in space and tab is the 
separation of the two objects; and ua and uab are unit 
vectors pointing along ta and tab, respectively. A new 
quantity, s ^ ta + tb has been introduced, in ad- 

dition. Inserting these expressions into Eq. (jB7p gives us 
SH^"^'' ^a, and the momentum of body A is then found by 
performing a surface integral over ^'s surface. The sur- 
face integrals are computed under the same assumptions 
as in Section IV Al namely the separation of the bodies 
is much larger than their radii, and each surface of inte- 
gration is a sphere immediately above a body's surface. 
When they are computed, one finds the same result as 
Eq. dUT]), 



SpA = mA5vA + \^^^SAy- riAB- (B9) 



One can find the momentum for body B by exchanging 
A and B. 

As a consistency check, we can evaluate the system's 
total momentum by doing a surface integral at infinity: 



(BIO) 



The quantity 5H^°''^^ is exactly the same as above, from 
which one can find 
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Wtot = ■m-Adv-'A + 2 



This, combined with the fact that 



as well as Eq. (jB9[) . gives 
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'field - q„2 

AB 



(B8d) as found in Section |V3 



{A ^ B) (Bll) 



5ptot = 5pA + 5pB + Spfi, 



iold- 



(B12) 



(A ^ B), (B13) 



[1] F. Pretorius, Phys. Rev. Lett. 95, 121101 (2005). 
[2] A. Ashtekar and B. Krishnan, Living Rev. Rel- 
ativity 7, 10. URL (cited on 5 February 2009): 



|http://www.liv ingreviews.org/lrr -2004-10| (2008) . 
[3] L.B. Szabados, Living Rev. Relativity 7, 4. 
(2004). URL (cited on February 21, 2009): 



15 



|http: / /www.livi ngreviews.org/lrr-2004-4 , 
[4] C.-M. Chen, J.M. Nester, and R.S. Tung, Phys. Rev. D 

72, 104020 (2005). 
[5] M. Campanelli, C. O. Lousto, Y. Zlochower, and D. Mer- 

ritt, Astrophys. J. Lett. 659, L5 (2007). 
[6] M. Campanelli, C. O. Lousto, Y. Zlochower, and D. Mer- 

ritt, Phys. Rev. Lett. 98, 231102 (2007). 
[7] J. A. Gonzalez, U. Sperhake, B. Bruegmann, M. Hannam 

and S. Husa, Phys. Rev. Lett. 98, 091101 (2007). 
[8] J. Healy et al., Phys. Rev. Lett. 102, 041101 (2009). 
[9] L. D. Landau and E.M. Lifshitz, Classical Theory of 

Fields (Addison Wesley, Redding Mass., 1962), Sec. 100. 
[10] C. W. Misner, K. S. Thorne and J. A. Wheeler, Gravita- 
tion (Freeman, San Francisco, 1973). 
[11] S. V. Babak and L. P. Grishchuk, Phys. Rev. D 61, 

024038 (2000). 

[12] F. Pretorius, in Relativistic Objects m Compact 
Binaries: From Birth to Coalescence, edited by 
Colpi et al.. Springer Verlag (see Sec. IV. C). 
|http : //xxx . lanl . gov/abs/0710 . 1338] 

[13] A. H. Shapiro, Vorticity, a film in a series pro- 
duced by the National Committee on Fluid 
Mechanics (1961); available for streaming at 
[http : / / web . mit . edu/hml/ncf mf . html | segment 4 



minutes 26 seconds into Part IL 
[14] K. S. Thorne and J. B. Hartle, Phys. Rev. D 31, 1815 
(1985). 

[15] J. D. Kaplan, D. A. Nichols and K. S. 

Thorne, Phys. Rev. D to be submitted, 

'http : //xxx . ianl . gov/abs/0808 72510 
[16] M. E. Rati and C. M. Will, Phys. Rev. D 62, 124015 

(2000). 

[17] R. D. Blandford and K. S. Thorne, Applications of Classi- 
cal Physics, textbook nearing completion and available at 
'http : / / www . pma . caltech.edu/Courses/phl36/yr2008 . 

[18] see e.g.. Appendix A of L.E. Kidder, Phys. Rev. D 52 
821 (1995). 

[19] T. Damour, P. Jaronowski and G. Schafer, Phys. Rev. D 

77, 064032 (2007). 
[20] K. S. Thorne, Rev. Mod. Phys 52, 299 (1980). 
[21] L. Blanchet, G. Faye, and B. Ponsot, Phys. Rev. D 58, 

124002 (1998). 

[22] H. Tagoshi, A. Ohashi, and B. J. Owen, Phys. Rev. D 

63, 044006 (2001). 
[23] G. Faye, L. Blanchet, and A. Buonanno, Phys. Rev. D 

74, 104033 (2006). 



